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We study the problem of body-force driven shear flows in a plane channel of width i 
with free-slip boundaries. A mini-max variational problem for upper bounds on the bulk 
time averaged energy dissipation rate e is derived from the incompressible Navier-Stokes 
equations with no secondary assumptions. This produces rigorous limits on the power 
consumption that are valid for laminar or turbulent solutions. The mini-max problem is 
solved exactly at high Reynolds numbers Re = Ut/v, where U is the rms velocity and 
u is the kinematic viscosity, yielding an explicit bound on the dimensionless asymptotic 
dissipation factor (3 = ei/U^ that depends only on the "shape" of the shearing body force. 
For a simple half-cosine force profile, for example, the high Reynolds number bound is 
P < TT'^/VSIG = .6715 . . .. We also report extensive direct numerical simulations for this 
particular force shape up to Re « 400; the observed dissipation rates are about a factor 
of three below the rigorous high-_Re bound. Interestingly, the high-i?e optimal solution 
of the variational problem bears some qualitative resemblence to the observed mean 
flow profiles in the simulations. These results extend and refine the recent analysis for 
body-forced turbulence in J. Fluid Mech. 467, 289-306 (2002). 



1. Introduction 

Bounds on the energy dissipation rate for statistically stationary flows belong to the 
small class of rigorous results for turbulence that can be derived directly from the incom- 
pressible Navier-Stokes equations without introducing any supplementary hypotheses or 
uncontrolled approximations. Quantitative approaches are mostly based on variational 
fori nulations as have been used in a variety of boundary-dri v en turbulent flows; see , 
e.g.,! Howard' ('1972|). IBuss3 fl978") , Doer ing fc Con stantin (1994") ,'Nicodemus et ai ('1998') 
iKerswell (,1998.) . More recentlv .Childress. Kerswell &: Gilbert. (2001) and .Doering fc Foiaa 
|)2002|) extended these analyses to body-forced flows in a fully periodic domain. The mo- 
tivation for such studies is to consider mathematically well-defined and tractable mod- 
els for (almost) homogeneous and (almost) locally isotropic stationary turbulence when 
boundaries are far away. 

Define the Reynolds number Re — U£/v, where U is the steady state rms veloc- 
ity, £ is the longest characteristic length scale in the body-force func tion and v is the 
kinematic viscosity. For the three-dimensional Navier-Stokes equations. IPoering & Foia3 
found that the energy dissipation rate per unit mass e satisfies 

e<cii^-^+C2— (1.1) 
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where the coefficients ci and C2 depend only on the functional shape of the body-force, 
and not on any other parameters or on any ratios involving the (say, rms) amplitude 
F of the force or the overall system size L — which could be arbitrarily larger than ^. 
(W e give a more p recise definition of the "shape" of the forcing function below, or else 
see lDoering fc Foiaa (,2002 ').') In terms of the dimensionless dissipation ratio /3 = e-^/C/"^, 
this result is /3 < Ci/Re + C2, an estimate in qualitative accord with theor etical, compu - 
tational, and experimental result for homogeneous isoptropic turbulence (Frisch 
ISreenivasan 1198 4), Sreeniyasan (1998)). The analysis in Childress, Kcrswell fc Gilberg 
l)200l|) focuses on dissipation estimates in terms of the true control parameter for such 
systems, the Grashof number Gr = Fl^ jv^ ^ but those results are less easily interpreted 
in terms of conventional ideas for homogeneous isotropi c turbulence. 

In this paper we refine and develop the approach in iDoering fc FoiasI ()2002|) for the 
particular example of fiow between free-slip boundary planes driven by a steady volume 
forcing density. Homogeneous shear turbulence is of interest in its own right, and by 
setting a uniform direction of the force we simplify some of the analysis allowing us to 
improve the bounds and make quantitative comparison with direct numerical simulations. 
The analysis produces rigorous limits that are approached within about a factor of three 
even at the moderate values of Re (up to 400) that we are able to reach computationally 
(in these runs the usual Taylor microscale Reynolds number R\ « 100). 

The rest of this paper is outlined as follows. In the next section we describe the model 
in detail and introduce some notation and definitions. In section 3 we derive the mini-max 
problem for the dissipation rate bounds. Elementary analysis quickly produces estimates 
which are then refined to the optimal bound (within this variational formulation) in the 
limit of Re — )■ 00. The final section 4 is a comparison of the results with the computational 
data and a brief discussion of the results. 



2. Preliminaries 

The three-dimensional dynamics of the fiow is governed by the Navier-Stokes equations 
for an incompressible Newtonian fluid, 

+ (u- V)u + Vp = i/V^u + f, V-u = 0, (2.1) 

at 

where the velocity field is u(x, i) — e^Ux -\- GyUy + e^Uz, the pressure field is p(x, t), and 
the kinematic viscosity is v. The flow is in the slab, [0, L^] x [0, £] x [0, L^] with periodic 
boundary conditions in the x (streamwise) and z (spanwise) directions. On the bottom 
and top surfaces at y = and y — £, we take the free-slip boundary conditions Uy =^ 
and dyUx = = dyUz- 

The steady body-force shearing the fluid is taken to be of the form 

f (x) = F^ (I) . (2.2) 

The length scale i is the longest length scale in the forcing function. The dimensionless 
shape function (p : [0, 1] — > 7?. satisfies homogeneous Neumann boundary conditions with 
zero mean: 

0'(O) = O = 0'(1), fmdri^Q. (2.3) 
Jo 

Technically we require that is a square integrable function, i.e., (p G L^[0,1], but 
in practice we are interested in even smoother functions whose Fourier transforms are 
effectively supported on low wavenumbers. The amplitude F is specified uniquely for a 
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given f when we fix the normalization of the shape function by 

1= / HvfdT]. (2.4) 
Jo 

We also introduce the dimcnsionlcss "potential" for the body-force shape function via 

f (x) = V X (I) e,] . (2.5) 

Then <I> € iJ^[0, 1], the space of functions with square integrable first derivatives, so 
^ = — and, without loss of generality because </> has zero mean, it satisfies homogeneous 
Dirichlet boundary conditions, $(0) = = $(1). 

At sufficiently high forcing amplitude, a finite perturbation causes transition to turbu- 
lence and the imposed driving sustains the turbulent state assuring statistical stationarity 
of the turbulent flow. In the following (•) denotes the space-time average. Using the root 
mean square value U = \J (u^) of the total velocity field including both a possible mean 
flow and turbulent fluctuations — and the length scale I in the force, the Reynolds number 
is 

i?e = ^ . (2.6) 

The energy dissipation per unit mass is e = z^(|Vup) and we deflne the dimensionless 
dissipation factor /3 via 

e = /?— . (2.7) 

Our aim is to derive bounds on /3 as a function of Re and as a functional of the shape </> 
of the driving force. 



3. Bounds for the energy dissipation 

The calculation of upper bounds on /3 proceeds in two steps. First is the derivation of 
a variational expression and second is the determination of rigorous estimates for it. 

3.1. The variational problem 

From the averaged power balance in the Navier-Stokes equations, the energy dissipation 
rate per unit mass e isf 

€=u(^\Vuf)=F{^u,). (3.1) 

Another expression for the forcing amplitude F can be obtained by projecting onto the 
momentum equation. Specifically, we project the streamwise component of the Navier- 
Stokes equations onto a mean zero multiplier function ip e i?^[0, 1] (a function whose 
second derivative is square integrable) satisfying homogeneous Neumann boundary con- 
ditions ip'{0) = = '0'(1)- The multipHer function ip must not be orthogonal to the shape 
function (j); we consider only {(f>ip) ^ 0. It is also convenient to introduce the derivative 
of the multiplier function, ^ = ip' £ H^[0, 1], satisfying homogeneous Dirichlet boundary 
conditions ^'(0) = = ^"(1). The inner product of (j) and V' is the inner product of $ and 
i.e., = {(l)iP) ^ 0. 

f Strictly speaking wc arc also assuming tliat the long time averages exist and that this 
relation is an equality for the solutions, rather than just an incquahty. That is, for weak so- 
lutions of the Navier-Stokes equations it is only known that t < F (4>Ux)- These mathematical 
technicalities do not alter the ultimate bounds that we will derive in this paper. 
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Take the inner product of the Navier-Stokes equation (|2.1(l with ip{y/£)ex, integrate 
over the volume utihzing appropriate integrations by parts, and take the long time average 
to obtain the relation 

-^iv^'w.u^^ = (3.2) 

This may be solved for the strength of the applied force F which when inserted into 13.1|l 
yields 

m ■ 

While the force amplitude F is not explicitly displayed in H3.3|l anymore, it is implicitely 
present through the constraint that the root mean square value of the velocity field is U . 
Dividing by ji, we produce an expression for the dimensionless dissipation factor (3, 

Changing now to normalized velocities ue^ + vBy + we^ = U~^(uxex + UyGy + u^e^), so 
that (v? + + w^) = 1, and dimensionless spatial coordinates £~^{xex + TJ^y + ze^), 
and using the potential $ and derivative multiplier ^' we have the identity 

The upper bound Pb on the dissipation factor is obtained by first maximizing the right 
hand side of H3.5() over all normalized, divergence- free vector fields satisfying the boundary 
conditions, and then minimizing over all multiplier functions ^ e [0, 1] satisfying 
homogeneous Dirichlet boundary conditions. Thus for any solution of the Navier-Stokes 
equations, (3 < j3h where the variational bound (Sb is the solution of the mini-max problem 

MRe) = nnnmax ^ \^^) ' (^.6) 

Note that while we explicitly display the Reynolds number dependence of it also 
depends the shape of the applied force — but not independently on the forcing amplitude 
F; the ratios in (|3.4|) . H3.5I) and H3.6I) are homogeneous in both (p and $. 

3.2. Evaluting bounds 

From (|3.6|l it follows immediately that /3f,(i?e) is bounded b y a function of the form 
ci + C2/Re for all Reynolds numbers, the analog of the result in lDoering fc Foias (,200^ . 
To see this, choose any convenient smooth multiplier function (e.g., $ and observe 
that elementary Cauchy-Schwarz and Holder estimates (recalling the unit normalization 
of u) give 

{<f'u) < (02)1/2 ^ ^^^^^ ^ 1 i^^^^i ^ ^ ^^,2^1/2 ^3^^^ 

^ y6[o,i] 

so that 

^>(«^)=^''^''':Tr'r' ''''''' ^''^'11':°^''^ (") 

This simple analysis produces explicit expressions for the coefficients c\ and C2 in a bound 
of the form 13 < ci/ Re + ci , displaying their functional dependence on the shape of the 



Energy dissipation in body-forced plane shear flow 



5 



driving force. In the following we will quantitatively and qualitatively improve this upper 
bound, computing the exact solution of the infinite Re limit of the mini-max problem. 

To further estimate and evaluate (3b{Re)^ note first that the boundary conditions to- 
gether with incompressibility imply that the y-component satisfies v{y) = where the 
overbar means horizontal and time average. We decompose the x-component into a hor- 
izontal mean flow u[y) and a fluctuating remainder u = u — u. Then the terms in the 
numerator of the ratio for /3h reduce: 



{(j)u) = {(jm) , (*{tt;) , 

Let = (m^). The normalization for the velocity field is 

so the terms in (|3.9|l may be estimated 



(3.9) 
(3.10) 



1 



\{^uv)\<^- sup \^{y)\{l-e), \{^'u)\<{-^'^f'^i. (3.11) 



ye[o,i] 



Hence for any choice of 5", 



max 



< max 



1/2 



1 



^ sup \^i,iy)\(i-e) + i-{^''y^'^ 



ye[o,i] 



Re 



(3.12) 

It is easy to find the maximizing value of ^. It is the solution of a quadratic equation 
in the interval [0, 1] for sufficiently high values of Re, or else it is = 1 if 



Re < 2- 



(3.13) 



supy6[04] |^'(y)r 

When = 1, the maximizing velocity field is a steady plane parallel flow, namely the 
Stokes flow for the given applied force. The right hand side of (|3.13l) is > 4, providing, 
by a somewhat round-about derivation, a lower bound for the smallest possible critical 
Reynolds number of absolute stability of the steady plane parallel flow that is uniform in 
the shape of the applied shearing force. For the purposes of the discussion here, however, 
we use the estimates above to bound (3i, as 

,1/2 



Pb < niin 



= mm ■ 



1 



max^(l-^'^)- sup |'I'(?;)| 



1 



o<C<i 
1 



sup |*(2/)| 



27 ye[o,i 



ye[OA] 
1 



max 
o<5<i^ Re ^ ' 



/2\l/2 



/2\l/2 



(3.14) 



This leads to improved estimates, in terms of variational problems for an optimal mul- 
tiplier ^E" and for C2 in a bound of the form /3b < ci/Re + C2. As we will now show, 
the variational expression above for C2 is sharp at high Reynolds numbers. That is, the 
Re —^ oo limit of the extremization problem for the optimal ^! and (3h can be solved 
exactly. 
Define 



/3(,(oo) = minmax ■ 
^ u 



First we will evaluate /3h(cxD), and then we will prove that 

limsup/36(i?e) < /3fc(oo). 

Re — ^oo 



(3.15) 
(3.16) 
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We accomplish this through a scries of two lemmas and a theorem. 

Lemma 1: If = e L2[0, 1] and * e iJi[0, 1] satisfies ^'(0) = = *(I), then 

max($'M) (^-uu) = ^(0^)^^^ sup |*(y)| (3.17) 
V27 2;e[o,i] 

where the velocity fields u are divergence- free and unit normalized in L^, periodic in 
downstream (a;) and spanwise (z) directions and free-slip in the normal (y) direction. 

Proof: We may take 5* 7^ 0. The calculation for (|3.14|l already established that the 
proposed answer is an upper bound to this variational problem, so all we must do is 
display a sequence of acceptable test fields u that approach the bound. Note that any 
nonvanishing ^ g H^[Q, 1] satisfying the homogeneous Dirichlet conditions is uniformly 
continuous and its extremum is realized at a (not necessarily unique) point j/m in the 
open interval (0,1). Consider the unit- normalized divergence-free vector field with 
components 

..=.g.(y)^/2sinfcz4--i=J|= 
Vk = ±gk{y)V2sinkz 

Wk — — g'k{v)V2 COS kz (3.18) 
k 

where gk{yf' is a smooth approximation of a (5-function with compact support centered 
on Dm and normalized according to 

1 /2 

The wavenumber A: 7^ is adjustable, and for each value of k we have {(t>Uk) = ((/>^) 
Now {"^UkVk) — {'^UkVk) and we may concentrate UkVk (y) = ±gk{y)'^ as tightly as desired 
around by taking k large, in which case {^UkVk) \^{ym)\ (gl) ■ Moreover, (g^) 
1/3 as fc — > 00. Hence there exists a sequence of test fields for which {(puk) {'^UkVk) 
approaches the upper bound. QED. 

Lemma 2: If $ e -^^[0, 1] and sign [$(y)] has a finite number of discontinuities, then 

mm — — — — — —. (3.20 

*6_f/i[04] ($^') (|$|) 

The minimizing function is 'l'm(y) ~ sign[$(y)] which is not in i7^[0, 1], but it is the 
pointwisc limit of a sequence of functions in H^. 

Proof: Note that 

sup \^{y)\m) (3.21) 
yG[o,i] 

so the proposed answer is a lower bound to the minimum, and the function "^miy) = 
sign [<I'(2/)] saturates this bound. Then it is straightforward to mollify the finite number 
of discontinuties in to produce a sequence of functions converging pointwise to 
\E'm,. Let the number of discontinuities of 5'm(2/) be N, located in order at ?/„. We may 
smooth by introducing a finite slope near each ?/„ to produce the regulated function 
^s{y) sketched in figure^ The mollified 'i'siy) is linear with slope zL6~^ inside all the 
intervals 26 around each i/„ and within 6 of the ends of the interval at y = and 1. 
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Figure 1. Boundary layer regularizatioii for the optimal multiplier 'i'm- (a) Sketch of ^miv) 
with a finite number of jumps in [0,1] ai y — 0,2/1,1/2, . . . , 1. (b) Sketch of the mollified 'i'siy) 
that converges pointwise to ^'m(j/) for 5-^0. 



Because <I> G ^^^^[0,1], $ = 0{y^\y — near the isolated zeros ?/„ where jumps. 
Thus (^^-a) = (|$|) [I -0(6^^^)]. Although = for each y G [0,1], 

^siy) ^ ^^miy) asS^O. QED. 

Theorem: If $ G ^^^[0, 1] such that sign [<i>(2/)] has a finite number of jump discontinu- 
ities, then 



/3f,(oo) 



and 



/27 (|<i>|) 

PbiRe) < /3fc(oo) + 0(i?e-3/4) 



(3.22) 



(3.23) 

as Re — + 00. 

Proof: Lemmas 1 and 2 establish the value of /^^(oo) in (|3.22|) . To establish (|3.23|) . recall 
from the proof of Lemma 2 that ($^5) = (|$|) [l - 0{S^/^)] and = 0{S-^). 

Using these facts together with Lemma 1, 



j3b{Re) = minmax 
^ u 



< min max 
~ (5 u 



< min max 
<5 u 



($*5) 



< min 

5 



1 (^^) 



1/2 



1 + 0{5- 



3/2^ 



^ /27 (|<i>|) 

Choosing (5 = ©(i^e^^/^) establishes the resuh. Qi^D 



(3.24) 



We make three short technical remarks here: 
(i) Although we only showed that hmsup^g^o^^ Pb{Re) < /?h(oo), it is natural to conjec- 
ture that at finite Re the optimal multiplier ^! actually looks like the mollified multipliers 
and that liuine^oo Pb{Re) = Pbioo). But this remains to be proven. 



8 C.R. Doering, B. Eckhardt and J. Schumacher 

(ii) The 0{Re^^/'^) rate of approach to the Re —f oo hmit in the theorem is not optimal 
for smoother shape functions. This is easy to see by repeating the proof of the theorem 
assuming, say, that </> e so that $ has a bounded derivative and behaves linearly 
(rather than as a square root) near its zeros. That generic hnear behavior leads to a 
faster 0{Re-^'^) rate. 

(iii) The hypothesis of a finite number of zeros in $ is probably not really necessary 
given $ G H^] we invoke it here for simplicity of the proofs only. In any case, for the 
applications we have in mind, 4> and <i> will actually be extremely smooth (composed, for 
example, of a finite number of Fourier components) so the theorem as stated and proved 
here serves our purposes. 

4. Comparison with numerical results and discussion 

Direct numerical simulations (DNS) in this geometry with these kinds of forces are 
possible in Fourier space thanks to the free-slip boundary conditions. For comp uta- 
tions we used the pscudospectral code developed in ^Schumac her fc Eckhardt (|2pOCf) and 
ISchumacher fc ICckhardt (2QQ1 ) with numerical resolution of 256 x 65 x 256 grid points. 
The steady volume forcing density, f(x), was chosen such that a laminar (and linearly 
stable!) shear flow profile uo(x) — — C/q coa{TTy / £)ex could be sustained. From the Navier- 
Stokes equations (|2.1|) it follows for this plane-parallel shear flow that 

f (x) = cos{ny/e)ex = F4> (|) (4.1) 

with shape <j){T]) = \/2cos7r77 and amplitude F = —vV n-K^ I . The flow can be consid- 
ered a Kolmogorov flow (see iBorue fc Orsza j lll996l) and lChildress. Ke rswell fc GilbertI 
l)200l|) ') with additional symmetry constraints in the normal (y) direction. The aspect 
ratio and scales for the calculations were L^j^ — L^/l — 2tt. The Grashof number for ab- 
solute (energy) stability of the steady plane-parallel flow with this force shape is Grc = 68 
where the Reynolds number is less than 7; the simulations were carried out well above 
this value, for Gr varying between 4900 and 59200. 

Mathematical results for the shape function 0(y) = — \/2cos7ry (equivalently <&(?/) = 
7r^^\/2sin7ry) are shown in Figure 2 along with the DNS results. In contrast to shear 
flows driven by rigid walls where the friction (dissipation) factor tends to decrease with 
increasing Reynolds number, here we observe a slight increase. The numerical values for 
P are about a factor 3 below the upper bound. This is a significantly better comparison 
of the data and the bounds than for turbulent Couette flow where the discrepancy is a 
factor 10 at Re « 10^. 

The mean profiles of the streamwise velocity component for two different Reynolds 
numbers are shown in the inset in Figure 2. It is interesting to note that even though the 
force shape is nonlinear across the layer, the mean profile is relatively linear with mean 
shear nearly constant outside boundary layers near the no-slip walls. The high Reynolds 
number limit of the optimal multiplier ipmiu) is piecewise linear with constant magnitude 
of its slope function; |'0m(y)| = |^'m(y)| — 1 away from the corners. We point out the 
similarity here with the observed mean profiles for the single example we have at hand. 

It will be very interesting to study the bounds (3h[Re) as well as the optimal multiplier 
functions at finite Reynolds numbers for a variety of force shape functions (/). This is 
because while the behavior of the boun d on (3 is similar i n structure to the observed 
experimental and computational values ijSreenivasa 1 lll984lll99a[) 1. it remains an open 
question how the high-i?e value of (3 depends on the details of the driving. There are 
some features we can anticipate right away, though. Assuming that the structure of 
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Figure 2. The dissipation factor {3 — et/U'^ as a function of the Reynolds number Re = U(./v 
for the force shape function (j){y/£) = —V2 cosTvy /£. The resuhs of the direct numerical sim- 
ulations are indicated by diamonds with error bars due to standard deviation zta/j where 
o"3 = /^(crf/e + 'i'^Urma /Urms)- The lower dotted line is the dissipation in the steady lami- 
nar flow which is a lower limit to the dissipation factor for any (statistically) stationary flow 
with this force. The three upper bounds, from top to bottom, are the estimate in 13.141 with 
the choice ^ = <1> (dash-dot line), the estimate in 13.1211 with the exact maximization over ^ 
followed by the choice ^' = $ (dashed line), and the optimal value /3b(oo) = 7r^/\/216 from the 
theorem with the rigorous 0{Re~'^^^) approach added on (solid). The optimal bound for the 
infinite Re limit is a 22% improvement below the infinite Re limit of the bound with 'I' = <1>. 
The mean flow profiles u^{y) for the simulations with the smallest (dashed) and largest (solid) 
Re are shown in the inset. 

^m{y) persists for large but finite Reynolds numbers, the high-i?e optimal multiplier is 
a simple but interesting nonlinear functional of the shape function of the driving force. 
While the plane-parallel Stokes flow profile Ustokes{y) is a linear functional of the shape 
function. 



the (infinite Re) optimal multiplier comes from a curiously similar — but highly nonlinear — 
formula: 



This expression for displays bifurcations as a functional of the shape function (ft. 
That is, for some shape functions, variations in (jj may result in no change at all in the 
associated high-i?e optimal multiplier, while at other configurations small changes in (p 
can produce large changes in (such as the number of "kinks" in the multiplier profile). 
Whether or not this kind of effect reflects any features of high Reynolds number mean 
profiles for shear turbulence driven by other shaped forces remains to be seen. 

To summarize, in this paper we have derived and analyzed a variational mini-max 
problem for upper bounds on the energy dissipation rate valid for both low and high 




(4.2) 




(4.3) 
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Reynolds number (including turbulent) body-forced shear flows. We find that the maxi- 
mizing flow fields are characterized by streamwise vorticies concentrated near the max- 
imal shear in an auxiliary "multiplier" p rofile, analogous to the "background" p rofile 
utilized in 'Doering & Constantin' ('l994(), iNicodemus et a\ l)l998|) . iKerswelll l)l998|) and 
IChildres's. Kerswell &: Gilbert [2001) . We solved the optimal high Re mini- max prob- 
lem exactly and compared the results with data from direct numerical simulations for a 
specific choice of forcing. We observed that the high Re bound is only about a factor of 
three above the data, and also that the high Re optimal multipier shares some qualitative 
features with the measured mean flow profiles. Future work in this area will include inves- 
tigations for other force shapes, as well as the improvement of rigorous bounds bv exact 
numer ical evaluation and/or by the inclusion of a balance parameter fsee lNicodemus et a\ 
in the variational problem. Finally, we remark that although we have carried out 
the bounding analysis with a steady driving for the flow, there is no obstruction to the 
inclusion of time-dependent forcing in the model. 

The numerical simulations were done on a Cray SVlex at the John von Neumann- 
Institut fiir Computing at the Forschungszentrum Jiilich; we are grateful for that support. 
This research was supported by the Deutsche Forschungsgemeinschaft and the US-NSF. 
Some of this work was completed while one of us (CRD) was resident at the Geophysical 
Fluid Dynamics Program at Woods Hole Oceanographic Institution during the summer 
of 2002, and helpful and stimulating discussions with N. Balmforth, P. Constantin, J. 
Keller and R. Kerswell are gratefully acknowledged. 
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